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High-density carbon nanotube networks (CNNs) continue to attract interest as active elements in 
nanoelectronic devices, nanoelectromechanical systems (NEMS) and multifunctional nanocompos- 
ites. The interplay between the network nanostructure and the its properties is crucial, yet current 
understanding remains limited to the passive response. Here, we employ a novel superstructure con- 
sisting of millimeter-long vertically aligned singe walled carbon nanotubes (SWCNTs) sandwiched 
between polydimethylsiloxane (PDMS) layers to quantify the effect of two classes of mechanical 
stimuli, film densification and stretching, on the electronic and thermal transport across the net- 
work. The network deforms easily with increase in electrical and thermal conductivities suggestive 
of floppy yet highly reconfigurable network. Insight from atomistically informed coarse-grained sim- 
ulations uncover an interplay between the extent of lateral assembly of the bundles, modulated by 
surface zipping/unzipping, and the elastic energy associated with the bent conformations of the nan- 
otubes/bundles. During densification, the network becomes highly interconnected yet we observe 
a modest increase in bundling primarily due to the reduced spacing between the SWCNTs. The 
stretching, on the other hand, is characterized by an initial debundling regime as the strain accom- 
modation occurs via unzipping of the branched interconnects, followed by rapid re-bundling as the 
strain transfers to the increasingly aligned bundles. In both cases, the increase in the electrical and 
thermal conductivity is primarily due to the increase in bundle size; the changes in network con- 
nectivity have a minor effect on the transport. Our results have broad implications for filamentous 
networks of inorganic nanoassemblies composed of interacting tubes, wires and ribbons/belts. 



The need for scalable synthesis of robust carbon- 
nanotube based devices has sparked the recent interest 
in high density carbon nanotube networks (CNNs). A 
diverse range of applications have been successfully re- 
alized and include flexible nanoelectronics^^^, chemical 
and mechanical sensing^ , heat transport^ multifunc- 
tional membrane^iitti^, supercapacitorJiS and nanocom- 
posites'^. While it is well-known that the properties of 
the CNNs degrade in comparison to those of the indi- 
vidual single- walled nanotubes (SWCNTs), their expres- 
sion is actively controlled by changes in the underlying 
filamentous nanostructure and its interplay with the net- 
work properties. Current understanding remains limited 
to the passive response. For example, both structural 
rigidity of and transport across these networks vary dra- 
matically with the degree of alignment and the density of 
junctions between the nanotubes^^"^-. Little is known re- 
garding the active interplay between the nanoscale struc- 
ture within these assemblies and their properties, yet it 
is of critical importance in almost all envisioned appli- 
cations. Of particular interest is the response to me- 
chanical stimuli and its effect on the electronic trans- 
port-'^^ The mechanics is expected to be non-linear, 
similar in principle to the mechanics of nanotube suspen- 
sions^" as well as filamentous assemblies of semi-flexible 
biopolymers such as actin, collagen and fibrirpitSSl 



deed, recent studies on severely deformed and physically 
interconnected, aligned nanotube networks have uncov- 
ered a thermally stable, viscoelastic responsd2ll_ Unlike 
these and the biopolymer networks, though, the as-grown 
and/or assembled CNNs that we study here are not phys- 
ically cross-linked and the SWCNTs can easily slide at 
the inter-tube junctions. In that sense, these as-grown 
CNNs represent an extreme limit wherein filaments of un- 
precedented rigidity interact via considerably weaker van 
Der Waals forces. In this article, we focus on the mechan- 
ics of aligned CNNs, assembled into a novel superstruc- 
ture, which allows us to directly quantify the interplay 
between mechanics and transport in these networks. We 
complement the experiments with coarse-grained com- 
putations that uncover the underlying nanoscale mecha- 
nisms. 



I. AS-SYNTHESIZED AND STRAINED SWCNT 
NETWORKS 

The architectural motif for our study is a sandwich 
structure consisting of millimeter-long vertically aligned 
SWCNTs sandwiched between two PDMS thin films, 
henceforth referred to as S/CNN/S. The superstructure 
is shown schematically in Fig. lA; the inset shows the 



densified 



PDMS laufif 
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FIG. 1. Nanostructural characterization of the CNNs with the S/CNN/S superstructures. (A, D, G) Schematic iUustration 
of the (A) as-fabricated (labeUed "pristine, here and elsewhere), (D) densified and (G) stretched network. Ti/Au metal pads 
are integrated into the PDMS layers which contact the nanotube terminations, as indicated. The imposed densification and 
stretching strains are approximately the same, e « 10%. (inset. A) Side view of the forest of vertically-aligned SWCNTs that 
are used to assemble the superstructure. (B-C) Low- and high-resolution SEM images of the aligned CNN within the pristine 
S/CNN/S superstructures. (B-C) Cross-sectional SEM images, (B) low-resolution and and (C) high-resolution, of the assembled 
superstructure showing the as-grown vertically aligned SWCNTs capped by PDMS layers. (E-F, H-I) Same as in (B) and (C), 
but for densified and stretched superstructures. The scale bars are 200 nm and 100 nm for the low- and high-resolution SEM 
images, respectively. 



side view of the vertically-aligned forest of SWCNTs syn- 
thesized using an ethanol based chemical vapor deposi- 
tion (CVD) technique. The entire assembly is capped 
by PDMS layers through a multi-step contact polymer 
transfer method. The interplay between mechanics and 
electronic transport requires characterization before and 
after deformation, and to this end, we have integrated 
nanometer thick metal contact pads onto the PDMS lay- 
ers that are in direct contact with the nanotube termina- 
tions; see Methods and Supplementary Information (SI) 
for additional details on the fabrication process. 

The as-grown SWCNTs are firmly embedded into the 
PDMS layers (SI, Fig. S3), a natural consequence of 
the strong PDMS- SWCNT interfacesP. The relatively 
long lengths of the SWCNTs is crucial as it ensures that 
the effects of the PDMS layers are minimized at equilib- 
rium and also during their active response. Additionally, 
each constituent nanotube in the vertically aligned net- 
work makes contact with both PDMS layers. The chief 
benefit of the fully percolated network is that the PDMS- 
SWCNT contact area remains constant during deforma- 
tion. While the absence of any interrupted nanotubes is 
a departure from the nanostructures that we expect in 
typical CNNs^^T their effect on the overall network con- 



nectivity is minimized since i) the density of these inter- 
ruptions scales inversely with the SWCNT length and is 
exceedingly small for the millimeter-sized nanotubes em- 
ployed here, and ii) the effect on the mechanics is limited 
to strains beyond which the individual nanotubes become 
extended; these are quite large since the as-grown verti- 
cally aligned SWCNTs are intrinsically wavy. 

Figures [Tj3 and[Tp show low- and high- resolution SEM 
images of the assembled CNNs within the superstruc- 
ture. Macroscopically, the SWCNTs are aligned, yet this 
is not strictly true at the nanoscale, evident from wavy 
nanostructure seen in the SEM images (Fig. [l]C). The 
filaments are largely bundles of nanotubes that are peri- 
odically rippled and also interconnected with an average 
inter-bundle spacing of r « 10 nm. The variations in the 
nanotube curvature and their lateral assembly that con- 
tributes to the quenched disorder is quite typical of the 
dense nanostructures associated with as-grown vertically 
aligned carbon nanotube forests and has its origins in 
the strongly confined and spatio-temporally non-uniform 
catalytic growth of the individual nanot ub e^SSHSSI 

We investigate the active response of the assembled 
S/CNN/Ss with respect to two classes of mechanical 
stimuli, densification and stretching, as shown schemat- 
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ically in Figs. [Tp and [Tp. For the former, the entire 
S/CNN/S is immersed in isopropyl alcohol (IPA) and 
then dried via evaporation (Methods). The eventual den- 
sification occurs in two stages: i) an initial surface tension 
driven elastocapillary aggregation during immersion as 
well as evaporation, and ii) further densification, mainly 
during evaporation, in which the bundles on average stick 
together due to the van der Waals forces once the lateral 
distance the between the bundles approaches the inter- 
graphitic spacin^23MI xhe net strain due to densifica- 
tion is e « 10% and the aligned CNN is transformed into 
rigid structure shown in Fig.jTJi;. The high-resolution im- 
age (Fig. [ip) shows that while the CNN is stiU textured, 
it is heavily interconnected with a slightly larger average 
bundle size. 

Our approach to impose stretching strains on the 
CNNs is fairly straightforward and illustrated in Fig. S2. 
The capping layers are pulled apart in a strain-controlled 
fashion using simple straining apparatus. The strain is 
applied to Si02 wafers that are retained on the cap- 
ping PDMS layers post-transfer for this specific reason. 
The CNNs are deformed to a maximum displacement of 
250 /im, or a uniaxial strain that is approximately the 
same as that during densification, e « 10%.- The CNN 
stretches easily, characteristic of a floppy network that 
deforms primarily via soft modes. The lateral compres- 
sion is negligible, i.e. the Poisson's ratio « 0. Un- 
like the viscoelastic response of physically interconnected 
CNNiP, there is no strain recovery following deforma- 
tion. The forces required to deform the network are 
quite small and we therefore do not extract the force- 
displacement curve and limit the characterization, both 
structural and the transport, to the maximum imposed 
strain, e « 10%. The low magnification image shows 
little change, yet the nanostructure is quite different 
(Fig. [Tj[). The network topology becomes increasingly 
aligned and there is an increase in bundle size and the 
network porosity. The frequency of the interconnections 
between the bundles is also reduced. 



II. TRANSPORT CHARACTERISTICS 

A. Electro-mechanically coupled response 

Two-terminal current- voltage {I—V) characteristics of 
the pristine and deformed superstructures are extracted 
in order to quantify the electro-mechanically coupled re- 
sponse. The curves for one of the assembled superstruc- 
tures is shown in Fig. [2]A.. The response is linear for 
the entire voltage regime (±0.5 V) and representative for 
each of the fabricated samples. The reciprocal of the 
slope is the total resistance, sum of the contributions 
from the CNNs - their intrinsic resistance, Rint - and 
that of the contacts, Rc- The mean resistance, R, aver- 
aged over seven different fabricated samples, is plotted 
in Fig. [2J3. Since the number of nanotubes that contact 
the metal pads is conserved during the two deformation 



modes, the effect of the mechanical stimuli on the contact 
resistance can be safely ignored. Then, the variation we 
see is due to the change in the intrinsic resistance of the 
CNNs. It is roughly the same for both modes and is sub- 
stantially reduced compared to the pristine CNNs. The 
strain factor, a quantitative measure of the coupled re- 
sponse, is {AR/R)/e « 6. Unlike the resistance change, 
the intrinsic strain factor of the CNNs will be higher as 
the extracted value is modified by the contact resistance. 




FIG. 2. Electrical and thermal characterization of the pris- 
tine (green), densified (red) and stretched (blue) S/CNN/S 
superstructures. (A) The measured current-voltage (J — V^) 
curves. The curves for strained superstructures are almost 
identical. (B) The mean resistance is averaged over seven fab- 
ricated superstructures (see SI). The errors bars are smaller 
than the solid circles and not shown. (C) Histograms for 
the thermal diffusivities a for the pristine and modified sam- 
ples, for three different fabricated S/CNN/S superstructures. 
(D) Average thermal conductivities of pristine, densified and 
stretched SWCNT networks. 



B. Thermo- mechanically coupled response 

The room temperature thermal characterization is 
based on thermal diffusivities extracted for pristine and 
modified CNNs. A steady-state temperature profile is 
established across the metal pads and directly yields the 
diffusivity a (see Methods). Figure shows the ex- 
tracted thermal diffusivities for three different fabricated 
S/CNN/S superstructures. In each case, the diffusivity 
decreases; a modest reduction for densified CNNs (34% 
on average) while an order of magnitude reduction for the 
stretched CNNs. The thermal conductivity, the prod- 
uct of the diffusivity and the volumetric heat capacity 
k = {pCp)a where p and Cp are the network density and 
specific heat respectively, is a more direct measure of the 
thermo-mechanically coupled response and is plotted in 
Fig. [2p for the three different CNNs. The conductivity 
of the pristine CNNs is well below those of individual 
nanotubes and their crystalline rope^^ even after we 
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account for the volume filling fraction of the network. 
The significant reduction is characteristic of the strong 
phonon-phonon Umklapp scattering that is highly sensi- 
tive to the nanostructure, in particular its network con- 
nectivitjEl. More importantly, the order of magnitude 
increase in the mean thermal conductivity suggests that 
the imposed deformation suppresses the scattering, ei- 
ther due to changes in network connectivity and/or the 
extent of bundling. The thermal factor associated with 
the thermo-mechanical coupling, {Ak/k)/e « 1. 



III. COARSE-GRAINED COMPUTATIONS 

The coupling between CNN mechanics and trans- 
port indicates an interplay between network density, the 
quenched disorder and extent of lateral assembly of the 
individual nanotubes, which we now fully explore us- 
ing experimental-scale computations. All-atom calcula- 
tions are clearly prohibitive due to the involved time- and 
length-scales. The inter-tube interactions are evidently 
central to accurately capturing the overall network re- 
sponse. To this end, we develop a coarse-grained ap- 
proach wherein i) the individual nanotubes are approxi- 
mated as linear chain of beads connected by springs that 
can bend, twist and stretch and ii) the van der Waal inter- 
tube interactions are extracted self-consistently using the 
universal graphitic potential integrated over the surface 
conformations of the individual nanotube^^SHMI ^j^g 
stiffness of the springs that connect the beads are tailored 
to the elastic properties of the individual nanotubes, ex- 
tracted from atomic-scale simulations and available as 
radius dependent empirical expression^sH Pqj. simplicity, 
we assume that the entire network is composed of nan- 
otubes of the same chirality and therefore radius, (8,8) 
SWCNTs with radius r = 0.543 nm. The inter-tube in- 
teraction is the sum of the interaction between beads on 
different nanotubes is modeled as a 6-12 Lennard- Jones 
pair potential. Classical many-body dynamics technique 
is employed to generate the conformations of the SWC- 
NTs within the initial equilibrated structure and to fol- 
low their evolution during the two classes of mechanical 
stimulP^. The resultant spring-bead approximation al- 
lows access to the structural and energetic evolution of 
a relatively large network consisting of millions of atoms 
and over several microseconds. 

Figures |3]^ and[3j3 shows the top and side views of the 
coarse-grained (bead) structure of the relaxed, percolated 
network that is part of a much larger CNN. The precur- 
sor as-constructed CNN is composed of 400 SWCNTs, 
each w 400 nm long and spaced 10 nm apart to match the 
inter-bundle spacing observed in the experiments. We fix 
the terminal ends of the nanotubes, a simple approxima- 
tion of the strong PDMS-SWCNT interfaces that stabi- 
lize the individual SWCNT ends within the S/CNN/S 
superstructure. The quenched disorder is generated by 
subjecting each nanotube to an initial compressive strain 
of sa 25%, a measure of the extent of disorder, and then 



dynamically relaxing the entire network to a local en- 
ergy minimum.'^ The SWCNTs immediately assemble 
into randomly branched and interconnected bundles com- 
posed of wavy and entangled conformations, quite like 
the network structure seen in the experiments. The lat- 
eral assembly is significant; the maximum bundle size 
of the network after relaxation is D^^ax = 53 nm. Fur- 
ther assembly is limited by the elastic energy cost associ- 
ated with bending, twisting and stretching the individual 
tubes; the computations allow us to directly quantify the 
energetic interplay between gain in the van Der Waal 
energy due to their self-assembly U^dw, and the elastic 
energy cost, Ue- 

Top and side-views of the coarse-grained nanostruc- 
ture within the stretched and densified networks (e = 
25%) are shown in Figs. [3]C-|3p and[3^-[3f^, respectively. 
The strain rate during the deformation is held constant, 
« 10^//is. The network is not allowed to compress later- 
ally as in the experiments (t^ ~ 0). The network poros- 
ity evolves as expected - decreases during stretching and 
increases during densification - yet the extent of the as- 
sembly visibly increases. We quantify the overall changes 
in the network by monitoring the two main parameters 
associated with its nanostructure, i) the maximum bun- 
dle size Dmax over multiple cross-sectional {XY) slices 
through the network (see Methods), and ii) the stored 
inter-tube interaction energy, U^dw i which is directly 
available. Their strain dependence is plotted in Fig. [Sp- 
[3)4. During densification, the bundle size increases mono- 
tonically and non-linearly and provides quantitative con- 
firmation of the enhanced bundling observed in the ex- 
perimental images (Fig. [l^). The increase in the magni- 
tude of the inter-tube energy (Fig.|3|4) is again consistent 
with lateral assembly of the bundles. We observe a simi- 
lar increase in the bundle size and |J7i,dVF| in the stretched 
CNNs but the temporal trends show a slight departure 
during early stages of deformation (e < 5%). Specifically, 
the network initially debundles; the average bundle size 
reduces by « 7% at an applied strain of e < 5% and 
iC^-udwl accordingly decreases. 

The network topology following the two deformation 
modes is qualitatively different, as seen in Figs. [3p-|3j^. 
During stretching, the entangled conformations of the 
individual nanotubes are ironed out during the initial 
debundling regime. The strains are accommodated by 
the interconnections that connect the branch points be- 
tween the entangled bundles. As a result, the bundles 
themselves are effectively shielded from the applied load 
and the interconnections are on average pulled apart that 
causes the bundles to dissemble. Past a critical point, the 
load transfers to the bundles and they start to assemble 
again, also evident from the trends in the inter-tube en- 
ergy (Fig. |3j4). More importantly, the bundles become 
increasingly aligned, which facilitates their assembly, and 
the network connectivity decreases. 

The insets in (Figs.[3pl and[3|l)l-[3]D3) highlight some 
of the nanoscale mechanisms that are involved in the 
network reconfiguration during stretching. Figure |3pl 
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FIG. 3. The strain dependence of the network nanostruc- 
ture in the computations. (A-E) Top (left) and side (right) 
views and of the coarse-grained network nanostructure ex- 
tracted from simulations performed at 273°K, for (A-B) the 
pristine, relaxed network, (C-D) network subject to e = 
25% uniaxial stretching strain, and (E-F) the densified net- 
work. The coloring scheme is based on the nanotube iden- 
tity of each discretized element, i.e. bead. The length of 
each SWCNT is 400 nm, the dimensions of the systems are 
200nmx200nmx400nm, and the average distance between 
SWCNTs is 10 nm. The circled regions highlight representa- 
tive mechanisms associated with the strain accommodation 
and their magnified view is also shown in panels CI, D1-D3 
and Fl. (G-H) Topology and energy change induced by defor- 
mation. (G) Evolution of the maximum bundle size Drnax and 
(H) the inter-tube energy Uvdw as function of strain during 
stretching and densification. 



and |3p2 show the zipping together of two smaller bun- 
dles into a larger bundle. The exact opposite, i.e. unzip- 
ping, can be seen in Fig. jSjDl wherein an interconnecting 
double strand (lower left) unzips from the surface of the 
bundle as it accommodates the imposed stretch. The 
combination of the two mechanisms, surface unzipping 
of a SWCNTs (colored green) and a double strand (col- 
ored yellow) and the zipping together of the two bundles, 
can be seen in Fig. |3jl)3. We do not observe substantial 



slip of the nanotubes and/or bundles since the network 
reconfiguration occurs mainly via soft modes. 

During densification, the increasingly bent conforma- 
tions of the SWCNTs and the bundles must also accom- 
modate the decrease in the average inter-bundle distance. 
At sufficiently large strains, the combination of lateral 
confinement and the branched network topology drives 
the lateral assembly. There is a concomitant increase in 
the alignment and network connectivity relative to the 
pristine network, visible from the network configuration 
shown in Figs. [3j5-[3|i^. The network reconfigures via sur- 
face zipping and unzipping at the bundles, similar to the 
nanoscale mechanisms in effect during stretching. The 
enhanced network density forces the nanotubes to pack 
more efficiently without dramatic increase in the elastic 
cost, and a dominant mechanism by which this is real- 
ized is the twisting of the bundles; one such representa- 
tive configuration is shown in Fig. [Sj^'l. In essence, the 
twist allows the network to store the additional increase 
in the energy stored by increasing the effective inter-tube 
interaction areaP^. 

Figure |4] shows the stress-strain curve during stretch- 
ing of the CNNs. The initial regime indicates a soft 
network with a very small yet positive Youngs modulus, 

« 20 ± lOMPa. This is characteristic of a fioppy net- 
work in that the connectivity is below the rigidity thresh- 
old for this class of systems. The strain accommodation 
is primarily through the changes in network connectiv- 
ity, lateral assembly and alignment, clearly visible in the 
network configuration corresponding to e = 35% (left 
inset). The modulus and the stress can also become neg- 
ative locally, an indication of transients associated with 
a self-healing network. At larger strains e > 35%, the 
distance between the fixed ends approaches the length of 
the SWCNTs and the network is no longer fioppy. The 
Young's modulus is of the order of ^ 50 GPa (indicated 
by the linear fit in Fig. |4]); the network configuration at 
e = 45% (right inset) reveals an almost fully formed bun- 
dle. 



IV. DISCUSSION AND CONCLUSIONS 

Our results challenge the traditional view of CNNs as 
a static network of carbon nanotubes and their bundles. 
Rather, a detailed analysis of the active mechanical re- 
sponse of the CNNs reveals a dynamic network that re- 
configures its nanostructure and the network topology via 
the extent of lateral assembly of the SWCNTs. This addi- 
tional degree of freedom modifies the network mechanics 
as it changes the distribution of the bundle diameters and 
therefore their stiffness. Clearly, existing network mod- 
els are ill-equipped to handle these strongly anisotropic 
and non-affine deformations modes and the relevant theo- 
retical modifications are addressed elsewhere [Wang and 
Upmanyu, in preparation]. The trends in the network 
reconfiguration during stretching and densification shed 
light on the electro- and thermo-mechanical couplings. 
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FIG. 4. Engineering stress-strain curve for the CNNs fol- 
lowing stretching. The inset shows the evolution of the total 
energy of the system. The snapshots of the network at large 
strains are also shown (e = 35% and e = 45%). 



which we now discuss in detail. The enhanced bundhng 
of the CNNs is the key to the observed increase in con- 
ductivities, in particular its effect on the electronic trans- 
port. Firstly, it is well-known that the electron disper- 
sion perpendicular to the bundle axis due to the inter- 
tube coupling decreases the separation of the conduction 
and valence bands around the Fermi level, in particu- 
lar the non-degenerate bands. The exact nature of the 
shift depends on the chirality of the nanotubes; it is as 
high as 0.25 eV for the bundles of (6,6) and (8,4) SWC- 
NT^'^^. Secondly, the doubly degenerate states also 
split within the bundle. As a classical example, strong 
band-splitting causes the (10, 0) nanotubes to become 
metallic upon bundling*''. The combination of these two 
effects reduces the band gap in semiconducting SWC- 
NTs. Of course, the reverse is also true for metallic nan- 
otubes, especially for armchair nanotubes, but the typ- 
ically higher ratio of semiconducting SWCNTs (the in- 
trinsic electrical heterogeneitjEl) in the as-grown forests 
leads to enhanced metallicity upon bundling. Lastly, the 
phonon transport that determines the thermal conduc- 
tivity also increases with bundling due to a non-linear 
increase in the number of transport channels along the 
bundle axi^. Interestingly, the non-monotonic variation 
in the extent of bundling observed during stretching of 
the CNNs (Fig. [sjGl) implies that the network initially 
undergoes a reduction in its conductivity during the de- 
bundling regime. Since this occurs at relatively small 
strains (e < 5%), it was not observed during the electri- 
cal and thermal characterization of the S/CNN/S super- 
structures. 

The individual SWCNTs are also strained during over- 
all deformation of the CNN, further modifying their 
band-gap. The strains, a combination of stretching, 
bending and twisting, can arise directly upon load trans- 
fer to the individual nanotubes, and indirectly due to the 



conformations of the individual SWCNTs at the branch 
points and within the bundles. For example, a net twist 
within a bundle occurs at the expense of bending and 
twisting of the individual tubes'^. Again, both chiral- 
ity and form of the strain modify the sign and extent of 
change in the band-gap, yet due to the higher density of 
semiconducting nanotubes in the heterogeneous forests, 
we expect the overall band-gap to become smaller. The 
conductivities are also modified by the network topology, 
in particular the extent of branching. In non-percolated 
networks, this is a major contribution as the metallic 
nanotubes are effectively shielded with increasing align- 
mentP^l. The effect is diminished here due to percola- 
tion, yet the branch points in a network serve as scat- 
tering sites for electronic and phonon transport. Once 
stretched well past the initial debundling regime, the 
topology evolves such that the branch-point density re- 
duces and further contributes to the increase in the ther- 
mal and electrical conductivity across the CNN. During 
densification, though, the network becomes increasingly 
branched as the nanotubes bend to accommodate the 
strains. Since the extent of bundling also increases, the 
overall increase in the conductivities is then primarily due 
to bundling and the debilitating effects due to scattering 
at the branch points are relatively less important. 

In summary, the combination the experimental char- 
acterization and experimental-scale coarse-grains compu- 
tations clearly show that the mechanics of carbon nan- 
otube networks is actively regulated by lateral assem- 
bly. Both densification and stretching result in enhanced 
bundling, although the network connectivity evolves dif- 
ferently. The increase in both electrical and thermal con- 
ductivities indicate that the bundling has a dominant ef- 
fect on the overall electron and phonon transport through 
the network. The extracted strain and thermal factors 
for the CNNs that we focus here are not exceptionally 
large yet have technological relevance because of the large 
strains that can be easily absorbed by the network. Al- 
though we have not explored the effect of average length 
of the SWCNTs, it is intuitive that increased lengths will 
enhance the connectivity as the individual SWNCTs can 
bend and twist more easily and the bundling will become 
more important. Our results also identify network scale 
features that can be engineered for controlled transport 
through these networks and have broad implications for 
a wide range of inorganic and organic filamentous net- 
works composed of assemblies of nanotubes, nanowires, 
and nanobelts. 



V. MATERIALS 
A. Synthesis of vertically aligned SWCNTs 

A multilayered substrate (Al/Si02/Si) is employed 
for the synthesis of the millimeter long SWCNTs. A 
20 nm thick Al buffer layer is first deposited on a 
SiO2(100 nm)/Si wafer. The surface is coated with a 
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0.7 nm thick Co catalyst layer using an e-beam evapo- 
rator and then placed inside a quartz tube. The tube is 
evacuated to « 15 mTorr and then exposed to a 5% H2- 
95% Ar mixture carrier gas at a pressure of 700 Torr and 
temperature 850° C. After the stabilization, high-purity 
anhydrous ethanol (99.95%) vapor is supplied from the 
bubbler as a carbon source which results in the growth 
of vertically-aligned SWCNTs. 



B. Assembly of the sandwiched hybrid 
superstructure, S/CNN/S 

A two step polymer casting method is use to synthesize 
the S/CNN/S superstructure. Each step is illustrated 
schematically in SI (Fig. SI) A bilayer metal contact pad 
(Ti(5nm )/Au(150nm)) is first deposited on the top- 
side of the as-grown vertically aligned SWCNTs using 
a sputter coater (Fig. Sla). A spin coater is used to 
coat the Si02 wafer with liquid phase PDMS (weight ra- 
tio of 10 : 1) and then w 80% cured. The top-side of 
the CNN with the metal pad is placed on the surface of 
the PDMS coating and then fully cured at 110°C for 90 
seconds (Fig. Sib). The Si02 wafers are detached from 
bottom (Fig. Sic and Sid) and the bilayer metal pad is 
deposited on the bottom-side of the CNN (Fig. Sle). It 
is finally capped with the second PDMS layer by pour- 
ing liquid PDMS and then curing at 110°C for two hours 
(Fig. Slf-h). 



C. Electrical and Thermal Chracterization 

A Janis ST-500 electric probe station connected to 
a Keithley 240 sourcemeter was employed for electri- 
cal characterization. The thermal conductivities of the 



three different SWCNTs networks were measured by us- 
ing standard laser flash thermal constant analyzer (Ul- 
VAC TC7000- HNC). 

D. Coarse-grained computations 

The radius dependent stiffness for SWCNTs are based 
on prior atomic simulatiorPSl. All stiffness used in our 
computations are for (8, 8) SWCNTs (stretching ks — 
626.46 eV/A, bending h = 9063.02 eV/A and torsional h 
— 6252 eV/A). The spring length connecting two neigh- 
bor beads has uniform length Al — 2nm and is small 
enough such that the orientation dependence of the in- 
teraction between beads can be ignored. The cohesive 
energy (—0.845 eV/nm) and the equilibrium inter-tube 
distance (1.4 nm) between parallel (8, 8) SWCNTs are 
readily available from the universal graphitic potential 
and serve as the two parameters for fitting the coarse- 
grained Le nnard - Jones pair potential, a ~ 1.534 nm and 
e = 0.81eM2il25l. The CNN is generated as described in 
the text and is stretched in a strain-controlled fashion 
by uniformly displacing the fixed ends of the SWCNTs 
(Z-direction). The Poisson's ratio in the experiments is 
negligible to simulate this, the cross-section area {X — Y 
plane) of the periodic computational cell is held constant. 
Densification is simulated by laterally compressing the 
computational cell while maintaining the distance be- 
tween the fixed ends of the SWCNTs. The time step 
for the dynamics is fixed. At — 5 is. To capture the 
changes in the network topology during physical load- 
ing, we track the bundle the maximum bundle size by 
sectioning the network into 1000 cross sections. Within 
each section, two SWCNTs are deemed to belong to one 
bundle if the inter-tube distance is smaller than the cutoff 
distance between the average first and second neighbors 
in an idealized hexagonally packed SWCNT bundle. 
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